arXiv:hep-ex/9808002vl 3 Aug 1998 


EUROPEAN ORGANIZATION EOR PARTICLE PHYSICS 


CERN-OPEN 98-020 
2 August, 1998 


Estimating the Parameters of Bose-Einstein 
Correlations from the Two-Particle Correlation 
Function in Multihadronic Final States 


A. De AngelisQ and L. VitaleQ 


Abstract 

To estimate the strength of the Bose-Einstein correlations and the radius of the 
hadronization region in multiparticle production, the two-particle correlation func¬ 
tions R for identical pairs is adjusted to a parametric function describing the en¬ 
hancement at small momentum differences. This is usually done by means of a 
binned uncorrelated least squares fit. This article demonstrates that this procedure 
underestimates the statistical errors. A recipe is given to construct from the data 
the covariance matrix between the bins of the histogram of the two-particle corre¬ 
lation function. 
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1 Introduction 


An enhancement in the production of pairs of pions of the same charge and similar mo¬ 
menta produced in high energy collisions was first observed in antiproton annihilations 
and attributed to Bose-Einstein statistics appropriate to identical pion pairs |^. 

Bose-Einstein Correlations (BEC) between pion pairs can be used to study the space- 
time structure of the hadronization source [Q. This has been done for hadron-hadron, 
heavy ion, muon-hadron and e+e“ collisions (see for example Ref. for reviews). 

The precise measurement of BEC parameters is especially important for LEP 2 physics. 
Interference due to Bose-Einstein correlations in hadronic decays of WW pairs has been 
discussed on a theoretical basis, in the framework of the measurement of the W mass 
1^: this interference could induce a systematic uncertainty on the W mass measurement 
in the 4-jet mode which is of the order of 40 MeV, i.e., comparable with the expected 
accuracy of the measurement. For a review of measurements of BEC in WW pairs, see 
for example p. 


2 Determination of BEC parameters 


To study the enhanced probability for the emission of two identical bosons, the correlation 
function R is used as a probe. For pairs of particles, it is defined as 


^(Pl,P2) 


P{Pl,V2) 
Po{PuP2) ’ 


( 1 ) 


where P{pi,p 2 ) is the two-particle probability density, subject to Bose-Einstein sym- 
metrization, pi is the four-momentum of particle i, and Po{pi,P 2 ) is a reference two- 
particle distribution which, ideally, resembles P{pi,p 2 ) in all respects, apart from the lack 
of Bose-Einstein symmetrization. 

If d{x) is the space-time distribution of the source, R{pi,p 2 ) takes the form 


R{Pi,P 2 ) = 1 + \G[d{x)]\‘^, 


where G[d{x)] = f dx is the Fourier transform of d(x). Thus, by studying 

the correlations between the momenta of pion pairs, one can determine the distribution 
of the points of origin of the pions. Experimentally, the effect is often described in terms 
of the Lorentz-invariant variable Q, defined by = (pi — p 2 )^ = M^Itttt) — 4m^, where 
M is the invariant mass of the two pions. The correlation function can then be written 
as 


R{Q) 


P{Q) 

PoiQY 


( 2 ) 


which is usually parametrized by the function 


R{Q) = n{i + . 


( 3 ) 


In the above equation, the pion source is spherically symmetric and gaussian, the pa¬ 
rameter r gives the RMS source radius, A is the strength of the correlation between the 
pions and N an overall normalization factor. The data from e’''e“ annihilations from PEP 
energies to LEP show values of r around 0.6 fm; the value of A strongly depends on the 
analysis technique. 
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It can be understood from what said above that the main problems in the study of 
BEC are given by a good choice of the reference sample, and by the definition of the 
normalization N. 

Most of the experimental studies are done inclusively, using all the observed charged 
tracks measured. Electrons, charged kaons and protons do not correlate with pious and 
clearly reduce the experimental correlation function. Another reduction of R{Q) is due 
to non-prompt pions, i.e., pions from decays of particles with lifetime larger than the 
hadronization scale of around 1 fm/c (like K°, A and b and c hadrons). These pions are 
not expected to correlate with those from the primary hadronization. 

3 Analysis 

The analysis performed in this paper was based on the JETSET simulation [§ . JETSET 
has shown to reproduce well the experimental inclusive distributions measured by LEP 

i- 

BEC were described by the Bose-Einstein simulation algorithm LUBOEI, fully inte¬ 
grated in the JETSET simulation. The values generated for the pion momenta are mod¬ 
ified by this algorithm, which reduces the differences for pairs of like-sign particles. This 
code has been shown [§ to reproduce well the two particle correlation functions measured 
in Z decays if Bose-Einstein correlations are switched on with a Gaussian parametrization 
for pions that are produced either promptly or as decay products of short-lived resonances 
(resonances with lifetime longer than the K*(890) lifetime were considered long-lived) and 
if the parameter values A = 1 and r = 0.5 fm are used as inputQ The value A=1 for direct 
pions corresponds to A ~ 0.35 for all pions or A ~ 0.25 for all particles [^. The htted 
value of the radius r depends on the choice of the reference sample. Using a Monte Carlo 
reference sample changes typically r from the input value of 0.5 fm to a 20% higher value. 

The use of simulated samples gives the possibility of knowing the normalization and 
of defining an unbiased reference sample simply by “switching off” the Bose-Einstein 
correlations. 


4 The naive approach to the fit gives wrong results 

We choose as a case study the hadronic decay of WW pairs with full Bose-Einstein effect 
and the BEC parameters set to A = 1 and r = 0.5 fm respectively. 400 samples of 2,000 
events each were simulated. For the reference sample Pq, 40,000 events were simulated 
without BEC, in such a way that the error on Pq gave a negligible contribution to the 
error on P. 

For each sample a histogram of the correlation function was built, using 40 bins of 50 
MeV each from 0 to 2 GeV. The normalization factor N was hxed by imposing that the 
average value of P between 1 GeV and 2 GeV was equal to unity. 

For each of the 400 samples we performed a fit to the form (^. The average values 
of A and r from the 400 fits were: 

A = 0.509 ±0.016, 

^The measured values of the parameters with a mixing reference sample for such “direct” pions in Z 
decays were A = 1.06 ± 0.17, r = 0.49 ± 0.05 fm |Q. 
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r = 0.569 ± 0.012 fm 

(the statistical error corresponds to the average of the errors). 

The “pull” of the htted values of A and r is shown in Figure ^ A Gaussian £t gives 
c’‘(a-<a>)/(ta ~ 1-35 with x^/DoF = 35/37, and (T{r-<r>)iar = 1-75 with x^/DoF = 30/38. 

The error on r is thus underestimated by a factor 1.75 while the error on A is un¬ 
derestimated by a factor 1.35. The fact that these factors are different indicates that 
the “naive” htting procedure described in this section could give possible biases in the 
determination of the parameters. 

When WW pairs in which one W decays hadronically and the other leptonically are 
simulated, the pulls are smaller (1.30 for r and 1.24 for A). This can be explained by the 
lower multiplicity, which gives smaller correlations. The case of the Z particle is similar 
to the semileptonic W. 


5 Improved technique for the evaluation 


The presence of bin-to-bin correlations in R is unavoidable: if there are, say, M positive 
tracks, the same positive track enters (M — 1) times in the two-particle density P. We 
build the generic nondiagonal term Cij [i ^ j) of the covariance matrix by assuming that 
it is given by the number of times that a track entering in the ith-bin enters also in the 
jth-bin. 

However, this is not the only statistical correlation effect, since the same track can also 
enter several times in the same bin. The latter effect can be accounted for by rescaling 
the error of each bin, cr(6j), to the expected relative error for the true independent counts. 
If we call Ti the number of extra counts due to tracks already present once in the bin i, 
then the number of independent counts di for the bin i is dj = 5* — r* and 


_ 1 
a / di 


(4) 


The diagonal terms Cu of the covariance matrix are assumed to be the squares of 

The correlation matrix pij = Cij/{aiaj) looks like in Figure |^a. In Figure the ratio 
between the errors on the bin contents estimated from Eq. (^ and the square root of the 
number of entries 6*. 

For each of the 400 samples we performed a x^ ht to the form (^, this time by using 
the covariance matrix described above. The average values of A and r from the 400 hts 
were: 


A = 0.499 ±0.021, 
r = 0.581 ± 0.017 fm 

(the statistical error corresponds to the average of the errors). 

The “pull” of the htted values of A and r is shown in Figure A Gaussian fit gives 
ct(a-<a>)/(7a = 1-02 with x^/DoF = 28/26, and cr(^r-<r>)/(7r = 1-09 with x^/DoF = 35/31. 

The error on r and A are thus correctly estimated (at better than 10%). A residual 
discrepancy can be attributed to the fact that, due to a Q-dependence of the fraction of 
particles which do non display BEG, the htting function (||) does not describe the data 
perfectly. 
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As an example, in Figure ^ the correlation functions R{Q) for like-sign pairs (closed 
circles) in one of the 400 simulated samples is shown together with the diagonal errors 
from the naive approach and from the approach proposed in this paper, and with the 
htting function corresponding to the average results of the £t in the two approaches. 


6 Conclusions 

Using a binned uncorrelated least squares £t to evaluate the parameters A and r of Bose- 
Einstein correlations underestimates the statistical errors and might bias the result. 

The error underestimate has been shown to be of the order of 30% in the case of single 
W and Z and of about 50% in the case of the hadronic decay of W pairs. A recipe is 
given to construct from the data a heuristic covariance matrix between the bins of the 
histogram of the two-particle correlation function. 
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Figure 1: (a) Pull function for the fitted parameter A using a binned uncorrelated least 

squares fit in the simulated samples. A Gaussian £t is superimposed as a solid line, (b) 
Same as (a) but for the parameter r. 
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Figure 2: (a) Correlation matrix (see text), (b) Ratio between the errors on the bin 

contents estimated with the procedure described in this paper and the errors computed 
from the naive approach. 
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Figure 3: (a) Pull function for the fitted parameter A in the simulated samples using the 

technique for the error estimate described in this paper. A Gaussian £t is superimposed 
as a solid line, (b) Same as (a) but for the parameter r. 



Q (GeV) 


Figure 4: The correlation function R{Q) for like-sign pairs (closed circles) in one of the 
simulated samples. The error bars shown are the square roots of the diagonal elements of 
the covariance matrix. The inner error bars correspond to the error computed from the 
naive approach. The dashed line shows the £t to Eq. (j^) where A and r are the averages 
of the htted values using the naive approach. The solid line corresponds to the approach 
proposed in this article. 
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